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Abstract. We compute in a relativistic way the time-of-arrival and the traversal 
time through a region of a free particle with spin i. We do this by applying the 
relativistic extension of the Event-Enhanced Quantum Theory which we have 
presented in a previous paper. We find a very good coincidence of the results of 
our formalism and the results obtained by using classical relativistic mechanics. 
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1. Introduction 

Blanchard and Jadczyk jl], ^, [| have proposed an extension of standard (non- 
relativistic) quantum mechanics called Event-Enhanced Quantum Theory (EEQT) 
which main idea is to view the total system as consisting of coupled classical and 
quantum part. The pure states of the quantum part are wave functions which 
are not directly observable, whereas the pure states of the classical part can be 
observed without disturbing them. Changes of the classical pure states are discrete 
and irreversible, they are called events. A review about applications of EEQT is for 
example [0. 

Blanchard and Jadczyk have also introduced a relativistic extension of EEQT |^| 
using the idea of a proper time and an indefinite scalar product. 

In a previous paper [|| , we have presented an alternative approach for a relativistic 
extension of EEQT. We summarize this approach in Section ||. 

The question when a particle arrives at a given point cannot be answered 
unambiguously in the standard formulation of quantum mechanics. Nevertheless, 
there exists a lot of approaches to answer the question of "time-of-arrival" . A review 
about "time-of-arrival" can be for example found in Muga, Sala, and Palao |7]] and an 
extensive review including a lot of references is Muga and Leavens B . 

In this paper, we consider a two dimensional spacetime. We examine the following 
operational definition of the time-of-arrival: The particle is prepared in a space-time 
point (0, xo) and moves freely (except of the influence exerted on it by the detector) in 
positive x direction. A detector is put at xd with xd > £o- It measures the time-of- 
arrival of the particle at xd- Because it is possible that the particle is never detected, 
the experiment or simulation should be stopped after a reasonable and finite period 
of time. 

This definition of time-of-arrival has been examined in the framework of non- 
relativistic EEQT [§. In Section | of this paper, we will compute the time-of-arrival 
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using our relativistic extension of EEQT and we will compare the results to those 
obtained by using classical relativistic mechanics of a point particle. 

The question how long a particle needs to traverse a given finite region cannot also 
be answered unambiguously in the standard formulation of quantum mechanics. This 
time is called traversal time. It is often examined if a potential is in the given region 
and the particle must tunnel through the region. Reviews about "traversal time" and 
"tunnelling time" including a lot of references are for example @ 0, @ 

"Tunnelling time" has been also examined using the non-relativistic EEQT 

§§!• 

In Section g of this paper, we examine the "free traversal time" for an always 
freely moving particle (again except of the influence exerted on the particle by the 
detectors). 

Using two detectors at rest one behind the other, we can measure the traversal 
time through the region located between the two detectors: The particle is again 
prepared in a space-time point (0, xq) and moves in positive x-direction. We put 
a detector D\ at x% with x\ > xq. This detector can detect the particle without 
destroying it. A second detector D2 is put at X2 with X2 > x%. It destroys the particle 
after detection. 

At the beginning of the measurement, both detectors D\ and D2 are active. If 
detector Di detects the particle (without destroying it), it turns itself off, but detector 
D 2 stays turned on. If detector D 2 detects the particle, the experiment is finished. 

Thus, the particle can be detected by detector D\ at a time t\ and then by 
detector D2 at a time t-2,. If this happens, the time difference t% — t\ is defined to be 
the "traversal time" . 

It is also possible that the particle is detected by D2 without prior detection by 
Di, but this situation should not contribute to traversal times. 

Moreover, it is possible that the particle is never detected or only detected one 
time by detector D\ . For this reason the experiment or simulation should be stopped 
after a reasonable and finite period of time. 

We will compute the traversal time using our relativistic extension of EEQT. 
Again we will compare the results to those obtained using classical relativistic 
mechanics of a point particle. 

2. A Relativistic Extension of EEQT 

We recall the extension of EEQT which we have proposed in || . It describes one spin 
i-particle in a relativistic way and in four dimensional spacetime. Here, we restrict 
ourself to consider two dimensional spacetime. 

As in EEQT, we postulate that the total system consists of a classical and a 
quantum part which are coupled. Because of that, at a given proper time r, the 
(pure) state of the total system is a pair {uj t , ^ T ). uj t is the state of the classical part 
and is the state of the quantum part. 

We assume that a (pure) state uj t of the classical part is a number: ui r G No = 
{0,1,2,...}. We also call a change of the classical (pure) state "event" as in non- 
relativistic EEQT. 

The (pure) states of the quantum part shall be (heuristically) solutions \P : 
K x K — > C 4 of the Dirac equation. Because we examine only free particles in this 
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paper, we use the free Dirac equation: 
d 

ihc-^-V(d,x) = H V(d,x) (1) 

with H a = — ic/ry^ 1 ^ + mc 2 7°. We denote the space of quantum states by 7i. 
A more precise definition of Ti. can be found in ||. We use the Dirac or standard 
representation of the 7-matrices. Let V = y 1 , a) : y Q ,y x ,a € R, \a\ < l} and 

<j\(u) = (y° + a ■ u, y 1 +u), Vu S M,VA = (y ,?/ 1 ,^) G P. We now introduce a 
positive-definite scalar product between two quantum states "J/ ,4, <J>_b 6 7Y: 

< * A |* fl > w := / j%d/„ (2) 

with A € f, = v I'^7 7 m 4'b and d/ M = (1, — a)du denotes the differential "surface" 
element of o~\. This scalar product is well defined because it is independent of A. 
This follows from Gauss theorem and the fact that d^j^g = 0. Moreover, one can 
show that this scalar product is covariant, its value being independent of the reference 
frame. 

We introduce the operators Uf ctotXo \ with do,Xo € R: 

(U {ctotXo) V)(x) := 9(ct , x + x) 

An interesting property of a quantum state is that it is uniquely given by its values 
on a space-like hyperplane o~\. Therefore, the operators U/ cto iXo \ are invertible. 
\& = U7^ q xo -,ip is the solution of the free Dirac equation ([!]) fulfilling the initial 
condition ^(do,x) = ip(x — Xq). We get 

*(ct,x) = (U[ c l o Xo) ^)(ct,x) = exp (~(t - t )H^ iP(x - x ) 

Now, we want to formulate an algorithm for modelling continuous relativistic 
measurements, indeed we will propose in the following an algorithm to simulate 
detections of the particle. In principle, we will do this by rewriting the algorithm 
of EEQT, replacing t with r and using our Hilbert space of "solutions." 

We denote the reference frame by K . The particle is prepared at proper time To 
in a point (do, Xo)- 

We consider n detectors with trajectories Zj(r), j = l..n. The trajectories start 
at proper time r = To from the backward light-cone of the space-time point of the 
"preparation event" : 

(d - z°(r )) 2 - (.t - z)(to)) 2 = 0, z°(r ) < d Q 

We allow detections which happen in the past of the preparation time. But we 
do not allow detections, if the detection space-time point is located in the backward 
light-cone of the space-time point of the preparation event. 

Each detector is characterized by operators Gj(r) : H — > H. Let Gt(t) be the 
adjoint operator. The total coupling between the quantum and the classical part is 
given by A(r) := E"=i ^(r^r). 

We define a detection algorithm in the following way: 

(i) The particle is prepared in a space-time point (do,Xo) at proper time r = tq. 
The quantum state is ^ To with H^toII-w =< ^toI^to >h= 1 an d the classical 
state is ui TQ — 0. 

(ii) Choose uniformly a random number r £ [0, 1]. 
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(iii) Propagate the quantum state forward in proper time by solving 

J^r = -iA(T)* T (3) 
until t = ti, where n is defined by 

1-11*^11^ = / dr < * T |A^ T > w =r 

Let lu t = lu to until t = tj , a detection happens at proper time r = Ti . 

(iv) We choose the detector k - which detects the particle - with probability 

Pk = ^\\G k (T 1 )* T1 f n 

with N = Y,UW G ^)*TA\n- 

(v) Let I be the detector which detects effectively the particle. The detection happens 
at the point zi(ti). The detection induces the following change of the states: 

The algorithm can start again perhaps with other detectors at position (ii). 

Because the scalar product is covariant, this algorithm is covariant. Moreover, 
its non-relativistic limit reduces to the algorithm of the non-relativistic EEQT. If we 
"charge conjugate" the initial state ^ To — > 4^ = C"f° T ^* o and the detector functions 

Gj(r) — > Cj° T G*(T) r y° T C + with C = i7 2 7°, then the algorithm will give the same 
detections as if we start with ty To and Gj(r) (if we choose the same random numbers). 
The quantum state in the "charge conjugated" world 4 , r 7 and the quantum state in 
the "normal" world are always connected by ^ = C7 W*. 

Note, that we have also formulated an algorithm for modelling ideal measurements 
of infinitesimal small duration in . It can be seen as playing the role of a relativistic, 
covariant reduction postulate. 



3. Initial quantum state 



The particle is prepared at proper time r = To in the space-time point (0,Xo) with a 
mean momentum p . We examine three different initial states of the particle in this 
paper. Remember that a initial states of the particle must be a solution of the Dirac 
equation ([!]). 

The first state corresponds to a state with only positive energies: 



dk- 



2E 



Po 



( E + m 



V k 



\ 



exp I \k{x — xq) — \Ect 



with rh= r -f,E = Vfc 2 + m 2 , Afc = 10 A \ 





F Ak (k) 



for \k\ < Afc 
otherwise 



and Np being a normalization factor so that H^o.pII^ = 1- This state describes an 
electron with charge — e. 
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The second one corresponds to a state with only negative energies: 



*o,Jv(ci, x) 
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exp ^— ifc(x — xo) + iEctj 
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I 2 



1. Remember that the above 



with Njy being a normalization factor so that ||\I r o,./v| 

algorithm is invariant under charge conjugation. If we consider the charge conjugate 
of the initial state and the detector functions, we get the same events. Because the 
"charge conjugated" world and the "normal world" should describe the same physical 
situation and because the charge conjugation of ^o.n describes a particle with charge 
+e in the "charge conjugated world" , we demand that the initial state ^o.w describes 
a positron with charge +e also in the "normal" world. 
As third initial state, we want to use a mixed state: 



0,PN 
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exp ^— ifc(a; — xq) + iEctj 



with 7] = 0.1 A. The constants are chosen in such a way that H'I'o.paHI^ = 1- ^ nc 
physical interpretation of the mixed state is the following: we assume that the particle 
(a single particle) can be in an "electron-state" (solution with positive energies) and 
in a "positron-state" (solution with negative energies), in analogy to the case, that a 
particle can be e.g. in a spin +i-state or in a spin — i-state. Superpositions as ^o.pn 
of the two states should be (in analogy to the spin-case) possible and allowed. 



4. Free Time-of- Arrival 



In this section, we apply the above algorithm to simulate the detection of the particle 
by one detector which is at rest. 

We want to compare the results to those which we obtain by using classical 
relativistic mechanics of a point-particle. 

Let us use the reference frame K in which the detector is at rest. In this 
reference frame, the particle is prepared at proper time tq = in the space-time 
point (0,x ) with a mean momentum p - The detector is put at sd, its trajectory is 
z(t) = (ct + x — xd, sd). The coupling operator should be given by 

G(r) = U^ T) g(x)U z[T) 

with g(x) characterizing the sensitivity of the detector: 

9{x) = y— ^— • F&xo(x) 
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The adjoint operator is G + (t) = U~^g + (x)U z (t)- 

Because it is possible that the particle is never detected, we stop the algorithm 
at r = tcut (with t C ut large). 

We want to recall that the algorithm is covariant. The choice of K as the 
reference frame is arbitrary. The algorithm can be applied in any reference frame, and 
it will result (if we choose the same random numbers) the same events in all reference 
frames. 

Using our algorithm, the probability that the detector detects the particle at all 
is given by 

/■TCUT 

Poc= dr < * r |A* r > H 

Jo 

The probability density for a "proper time-of-arrival" at the detector is given by 
(r < tout) 

P(t) = ^-< * T |A* r > H 

It is zero for r < and t > tcut- 

Using this probability density for "proper time-of-arrival" , we can calculate the 
probability density and the expectation value for the time-of-arrival in an arbitrary 
reference frame. 

Let us first look at the detector's rest-frame K . If a detection happens at proper 
time t, then it happens in space-time point z(t) = (cr + x — Xd,Xd)- This implies 
a time-of-arrival of t = t — XD ~ X ° . So we get the following probability density for the 
time-of-arrival in the detector's rest-frame Kq: 

x D -x Q 



Qo(t)=p{t + 

The expectation value (or mean time-of-arrival) is 

T a , = J dtt Qo (t) = y dr (r - p( T ) = / drr P (r) - 

Now, we want to calculate these values in a reference frame K v which moves with 
velocity v with respect to the detector's rest- frame Kq. The Poincare-transformation 
Ko — ► K v has the following form: 

x = 



i -4 v — 1 

The detector trajectory in K v is 

Z ( T ) = I 1 _ j ■ (ct + x - x D - v/c ■ x D , -vt - v/c ■ x a + v/c ■ X D + X D ) 

So the normalized probability density for the time-of-arrival in the reference frame K v 
is given by 




v 2 ~ xd — xq v \ 
Qv(t) = V 1 - — •Ph/l--ji + + I 

and the expectation value (or mean time-of-arrival) in K v is given by 

T a , v = ( dttg v (t) = r i [t Oi0 - ^x D ] (4) 

J v 1 - ^ 
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4-1- Numerical Approach 

We use the reference frame K to compute p(r). Therefore, we define 
fl(r, x) := {U z ^ T )^ T )(x) = * T (cr + x - x D ,x D + x) 
If ty T is a solution of ([]]) and (0), then we get 

ifi— ft(r, x) 

■ t 9 x / * > \ . t /9* r V 

= ifa ^pj W T (cr + x - x D ,x D +x) + in I ) ( cr + x - x D , x_d + x) 

= i/oO(r,a ; )-i^ + (x) 3 (x)fi(r, a; ) (5) 

We have to solve this equation with the initial condition O(0, x) = ^q(xq — Xd, Xd+x). 
Using < ty T \A$F T >-h = / dx fl + (t, x)g + (x)g(x)Q(T, x) , we can calculate and p(r) 
if we know CI(t,x). Using p(r), we get Qo(t) and T aj o- 

The equation (||) with the initial condition O(0, a;) = \Po (^o — %d , xd+x) is solved 
numerically. The proper time dynamics of Q is approximated by 

Q(r + Ar) « exp ( — ^-\g+( x )g{x) ] exp ( -Ar^ ) exp f ] fi(r) 



We now discretize the proper time and the space with steps Axb = cAtb = 0.0004 A. 
Then, the first and the last operator can be computed directly. The second operator 
is discretized by using the method of Wessels, Caspers, and Wiegcl JlT] . The 
boundary conditions are walls at x — — 6 A and at x — 4 A in such a way that 
n(r,-6A) = 0(r,4A) = for all r. We set: t CU t = 13.0A/c(p < 0.5 mc), 
tcut = 7.0A/c(0.5mc < pa < 0.75 mc), tqut = 5.0 A/c(0.75mc < pa < l.Omc), 
tcut = 4.5A/c(1.0 mc < po). We do the simulations again with other time and space 
steps Ax a = cAta = 0.0006 A. So the error in the expectation value T Qj o can be 
approximated by 

error(T afi ) = ±^-^^— \T afi {Ax B ) - T afl (Ax A )\ (6) 

4-2. Results 

We set x = -1 A, x D = OA, Ax D = 0.01 A and W D = 1 x lO^mc 2 . 

Figure |l| shows the corresponding expectation values of the time-of- arrival T Q) o 
in the detector's rest-frame K$ for different momenta po and for the three different 
initial states. The error bars are calculated using (Q). 

In addition, Figure |l] shows the arrival-times calculated by using the classical 
relativistic mechanics of a point-particle: 

t a ,RM 

We see that the expectation values are nearly independent of the initial state 
^o,pn, ^q.p or \&o,iv- Furthermore, it exists a good agreement between the values we 
computed and the results obtained by using classical relativistic mechanics. 

Only for very high momenta, the expectation values of the simulation with ^o.pjv 
are a bit smaller than the times from classical mechanics and those obtained by the 
simulations with other initial states. 
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Figure 1. Mean time-of-arrival T a ,o versus particle momentum po in the 
detector's rest-frame Kq, relativistic simulation with detector parameters Ai^ = 
0.01 A, Wd = 1 x 10~ 5 m<? started with different initial states : ^o,P (boxes with 
error bars), jv (triangles with error bars), PN (circles with error bars), other 
parameters see text; classical relativistic mechanics t a ,HM (dotted line); the figure 
inside is a zoom of the right lower area of the figure outside 



The reason can be seen in Figure which shows probability densities in the 
detector's rest-frame Kq. For very high momenta and if we start with ^o,pn, we find a 
small probability for negative times-of-arrival. This fact explains why the expectation 
values of the simulation with ^o.pjv are smaller than the results of classical mechanics 
and those of the other simulations. 

We also see that the probability densities are (nearly) the same if we start with 
^0,P or * ,JV- 

The expectation values T a ^ v in different reference frames are connected by (Q). 

We get from classical relativistic mechanics: the time-of-arrival t a ^M in the 
reference frame K v is connected to the result t a ,RM in the reference frame Kq in 
the same manner: 

v 

ta,RM ~XD 

. C l J 

V 1 c 2 

So we have a good agreement between the simulated expectation values and the 
results deduced from classical relativistic mechanics in all reference frames! 

Another important and interesting question is how the expectation values depend 
on the parameters of the detector. The initial state is now the function Vl/o.p with 
positive energies. We examine those particle momenta which are also examined in 
Figure |l|. We compute the probability densities and the expectation values for four 
different pairs of detector parameters. 

First, we examine the case of a "higher" detector with Wd = 1-0 mc 2 , but its 
width is still Axd — 0.01 A. We find that the expectation values and the normalized 
probability densities do not change by using a "higher" detector for all examined 
particle momenta. The detection probability increases with increasing detector 
height Wd- 

Next, we examine the case of a wider detector with Axd = 0.4 A and Wd = 



a.RM 
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Figure 2. Probability densities qq for the time-of-arrival in the detector's rest- 
frame Ko, detector parameters: Axp> = 0.01 A, Wo = 1 X 10 — 5 mc 2 , initial state: 
\l/o p (small solid line), n (big dotted line), ipg PJV (dashed line), particle 
momentum po; the vertical solid line indicates the arrival time deduced from 
classical relativistic mechanics 



1 x 10~ 5 mc 2 . The expectation values do not change. The normalized probability 
density do not change in a significant way, it only becomes a bit wider. Again the 
detection probability Pqo increases with increasing detector width Axd- 

The results only change if we use a very wide and height detector with Axd = 
0.4 A and W D = 1.0 mc . The expectation values and the normalized probability 
densities are then shifted to earlier times. 

In other words, the simulations show that it exists a wide range of detector 
parameters for which the results do not change significantly. 

5. Free Traversal Time 

In this section, we simulate the traversal time measurement described in Section [j] by 
applying our algorithm being described in Section ^. 

We use the detectors' rest-frame K . The particle is prepared at proper time 
r = in (0,xq). The first detector D\ is at rest at position x\. Its trajectory is 
z\{t) = (ct + xo — Xi,Xi). The second detector D 2 is at rest at position x 2 - Its 
trajectory is Z2(t) = (ct + xq — #2> Xi)- The coupling operators of detector Dj should 
be given by 

Gj{r)=U-\ T)9j {x)U Zj{T) j = 1,2 
with g(x) characterizing the sensitivity of the detector Dy. 

9j(x) = J^jr 1 -F^ix) 
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Because it is possible that the particle is not detected two times, we stop the 
algorithm at r = tqut (with tqut large). 

Let be the initial state and ^ T the solution of (|l|) and (|3|). Then, the 
probability that the particle is detected by D\ at all is 

rTCUT 



Poo,i = / C T dr < ¥ t |G+(t)Gi(t)¥ t >h 
Jo 



The probability density that the particle is detected by D x is given by (r < tcvt) 

Pi(r) = -^—< * T |Gfi-(r)Gi(r)* T >n 

If a detection by detector D\ happens at T\ , the quantum state after the detection is 
given by: 

(ri) >= gx(n)^ 

Let $r Tl ' ) be the solution of (^) with initial state (^). We get the following conditional 
probability that the particle is detected a second time by D 2 if it is detected by Di 
at Ti. 

P£ l] = r UT dr 2 <$^)|G+(r 2 )G 2 (r 2 )^) > n 

The probability density for a second detection at proper time t 2 by detector D 2 after 
a detection of detector D\ at proper time n is given by 

P 2 n) (r 2 ) = -^-y < $^|G+(t 2 )G 2 (t 2 )<I>£> > w 

Finally, the probability density for a first detection by D\ at n and a second detection 
by D 2 at r 2 is 

, , P 2 Tl) (T 2 )-Pic Tl) - Pl (T 1 )-P 00 , 1 

Pl2(Tl ' T2) = rc/T, p (r) — 7—^ — 

< $(; i )|G+(r 2 )G 2 (r 2 )<i>(.; 1 ) > w • < * Tl |Gt (n)Gi(n)* n > w 
1 I for < ti < tcut and t\ < r 2 < tcut 



otherwise 

with Poo,i 2 being the probability that the particle is detected two times: 

Pcv.ii = Poo.i ■ I drPW •pi(r) 
Jo 

= f CL dn / Cf/ dr 2 < $^)|G+(r 2 )G 2 (r 2 )$(.; 1 ) > H • < *n|Gf(ri)t3 ! i('n)*r 1 >W 

JO Jti 

Note, that this probability density is independent of the reference frame in which 
the algorithm is applied. 

We now calculate traversal times in different reference frames. In contrast to pi 2 , 
the probability density for traversal time depends on the reference frame. 

We start with the detectors' rest-frame Kq. If the first detection of D\ happens 
at proper time n, then it happens at space-time point zi(ri) = (ctx + xq — x%,xi). If 
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the second detection of D 2 happens at proper time t 2 , then it happens at space-time 

z 2 (r 2 ) = (cr 2 + Xq — x 2 ,x 2 ). The resulting traversal time is therefore 

xo - x 2 Xq — xi x 2 - Xi 

t = T2-\ n = T2 — T\ 

C C C 

So the normalized probability density for the traversal time in the detectors' rest-frame 
Kq is given by 

/ \ f ( X2 — Xi 
po{t) = / dr P12 [r,t-\ hr 

The expectation value of the traversal time (or mean traversal time) in Kq is 

T t , = J dtt- J dr Pl2 (r,t+ X2 ~ Xl +r 

Now, we want to calculate these values in the reference frame K v (the reference 
frame which moves with velocity v with respect to the detectors' rest-frame Kq). 
The detector trajectories in K v are 

_ . . 1 / v v \ 

Zi{t) = — • CT + XQ — X\ Xi, —VT (X — Xi) + Xx ) 

/-, v 2 V c c / 

1 / WW 
-£2, -UT ( 

c c 

If the first detection of Dj happens at T\ and the second detection of D 2 happens at 
r 2 , then it results a traversal time of 

r 1 ( x 2 -xi v 

t = , • r 2 - n ? (x 2 - £Ci) 

c <r 



Z 2 (T) = ■ \CT + X Q -X 2 X 2 , -VT [Xq - X 2 ) + X 2 I 

t/l-£ c c / 



1- s v 

So the normalized probability density for the traversal time in the reference frame K v 
is given by 



Pv (t) = ^/l-^ J drp 12 \r^l-—t+ — + -( X2 - X1 ) +T I 

The expectation value of the traversal time (or mean traversal time) in K v is 



Xp" ~ X2 — X\ V 

T t , v = I dtt ■ \/l - "— ■ I drpu [ t, \/l - -^t H 1- — (x 2 - xi) + r 

Tt.o-^a-^i) ( 8 ) 



1 - 



5.1. Numerical Approach 

We use the reference frame Kq. For computation of the algorithm until the first 
detection, we define 

fl^(r, x) := (C7'( CT+a . _ a . li o)* T )(a;) = * T (cr + x - xi,x) 

ty T should be a solution of ([!]) and (||), so we get 

d ft 
ih—n A (T,x) = H n A - i-gf(x - xi)gi(x - Xi)Cl A 

-i^Tg 2 (x - x 2 )+<? 2 (x - x^T-^a (9) 
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with T = U icT+Xo _ Xu0) U~l +XQ _ X2 fi) = exp (-(x 2 - x{]-^H Q ). A solution of this 
equation has to be found satisfying the initial condition f2^(0,x) — ^o(xq — x\,x). 
The equation (^) is solved numerical with the proper time dynamics approximated by 

^(t + At) ps 

/ Ar i 2 Ar 1 , A ( At i ( . Q x d 

° XP ft mC 7 " — 2 51 ( " " " " l} J 6XP {~2"R \~ lhC1 7 to 

Texp \-AT~g£(x- x 2 )g 2 (x- x 2 U T^ 1 exp (-^^ ( -tticj ^ 1 



2 h\ 11 dx 

f At i 20 At 1 \ / >\ n , > 

exp I — 2-^ mc 7 ~ — 2-9i (' T ~ ^iMO ~ ^l) J "a(t) 

with T re J] exp (-At^Hq) . 

We discretize the proper time and the space with steps Ar and Ax (cAr = Ax). 
The boundary conditions are walls at x = — 8 A and x = 8 A in such a way that 
fU(T,-8A) = 04(V,8A) = for all r. 

All operators (including T) can be evaluated directly or are approximated by 
using the method of Wessels, Caspers, and Wiegel |l7| or by using Wendroff 's formula 
(see e.g. @). 

For simulating the second part of the algorithm (after a first detection by detector 
Di at proper time ti), we define 

( T > X ) : = ( U (cr+xo-X2.,0)^r)(x) = * T {cT + X - X 2 , x) 

with ty T being a solution of ([j]) and (||). We get 

\K^P$\t,x) =H n£\T,x)~£g+(x-x 2 )g 2 (x-x 2 )n£\T,x) (10) 
We must solve this equation with the initial condition 

y / dxQ\(Ti,x)g^(x - x 1 )g 1 (x - xi)CIa{ti,x) 

The equation ( |l0| ) with the initial condition ([ll]) can be solved approximately in 
analogy to Section 4.1. 

Knowing £Ia(t,x) and f2jj (r, x), we can calculate Poo, 12 and P\ 2 (ti,t 2 ). So we 
can calculate po(*) and Pt,o- 

We do the computation with proper time and space step cAr# = Ax_b = 0.0006. 
The value of tcut depends on the particle momentum: tqut — 31.5A/c(po = 
0.25 mc), tqut — 17.5A/c(po = 0.5 mc), tcut = 13.5A/c(po = 0.75 mc), tqut = 
11.5 A/c(1.0mc < po < 1.5 mc), r C)7T = 10.5 A/c (1.5 mc < p )- 

Moreover, we do the computation with proper time and space step cAta = Ax^ = 
0.001. So the error in the expectation value T tj o can be approximated by 

error{T tfi ) = ± Ax * |T t , (Ax B ) - T tfi {Ax A )\ (12) 

Ax a — Axb 

The error in the probability Poo ,12 is approximated by a similar formula: 

erroriP^) = ± A |P 00 , 12 (Ax B ) - P 00 . 12 (Ax j4 )| (13) 

Ax a — Axb 
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Figure 3. Mean traversal time Ttfi versus particle momentum po in the detectors' 
rest-frame Kq, starting with different initial states : p (boxes with error bars), 
*5>0 N (triangles with error bars), V&o.PJV (circles with error bars), other parameters 
see text; results from classical relativistic mechanics t t (dotted line) 




5.2. Results 

We perform the simulation with different initial states and different particle momenta 
Pq. We set xo = —1.5 A. The detector parameters are x\ = OA, Axi = 0.5 A, 
Wx = 1 x 10~ 3 mc 2 and x 2 = 1.26 A, Ax 2 = 0.02 A, W 2 = 1 x 10~ 3 mc 2 . 

Now, we are interested in the traversal time in the detectors' rest- frame K . 
Figure ^ shows the expectation values for traversal time with different initial states 
and different particle momenta po- The errors calculated by (|l^) are also plotted. 

A first result is that we see nearly no dependence on the initial state. 

In addition, the times which one obtains by using classical relativistic mechanics 
of a point-particle are plotted: 

x 2 - Xi 

tt.RM = 

C 

There is a good agreement between the simulated results and those obtained by 
using classical relativistic mechanics. This agreement becomes more accurate with 
increasing particle momentum po . 

Figure || shows the probability densities po for traversal time in the detectors' 
rest-frame Kq with different initial states. The probability densities have a peak at 
the classical expected traversal time. Again, we see that there is nearly no difference 
if we start with the initial state ^o.p or ^o.w- There are only small differences with 
the results obtained with the initial state $ ,pjv- 

Next, we examine the situation in a moving reference frame K v . It moves with 
velocity v relative to Kq. We calculate the traversal time in the framework of classical 
relativistic mechanics: 

1 r v , ' 

U,RM = — , ' H,RM 7j\X 2 — X\) 

\ v 2 L c z J 

r 
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Figure 4. Probability densities po for traversal time in the detectors' rest-frame 
Ko, initial state: ^o,P (small solid line), ^o.iv (big dotted line), "fr^piv (dashed 
line), particle momentum po; the vertical solid line indicates the traversal time 
given by classical relativistic mechanics 



We get the same correlation between tt,KM and tt,RM m classical relativistic mechanics 
than the correlation between T t ,o and T t<v in our formalism (see (||)). We get the same 
result than in the time-of-arrival case: the good agreement between our results and 
those obtained in classical relativistic mechanics exists in all reference frames. 

Now, we want to examine how the results depend on the detector parameters. 
The particle momentum is fixed at p = 0.75 mc and the initial quantum state is ^o,p- 

We start by varying the parameters of the first detector D\. The parameters of 
the second detector Di are fixed at Ax2 = 0.02 A and W2 = 1 X 10~ 3 mc 2 . 

First, we compute Poo, 12 and the expectation value T tj o in Kq for different detector 
widths Axi while keeping W\ — 1 X 10~ 3 mc 2 fixed (see Figure |B|a). We find out that 
it exists a range of detector width (0.3 A < Ax\ < 1.0 A) for which the expectation 
value T t> o does not change in a significant way. But the probability for two detections 
-Poo, 12 increases with increasing detector width Ax±. 

In the range 0.3 A < Axi < 1.0 A the forms of the probability densities po do 
not differ in a significant way. The peaks only become wider with increasing detector 
width Ax\. If the detector width is very small (Axi — 0.02 A), the wave function 
changes strongly through the detection by D\ and so we get a qualitatively different 
probability density pg. 

Now we fix Axi — 0.5 A and vary W\ (see Figure |^b). In the case of weakly 
intrusive detectors W\ < 5 x 10~ 3 mc 2 , the expectation values T t> o do not differ in a 
significant way. For higher detectors, the expectation values T tj o increase a bit with 
increasing detector height W\. The probability Poo, 12 increases with increasing Wi, a 
fact one expects intuitively. 

With increasing detector height W%, the peak of the probability densities po is 
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Figure 5. Mean traversal time Ttfi (circles with error bars, left axis) and 
probability Poo, 12 (boxes with error bars connected with a solid line, right 
axis); initial state fop with po = 0.75mc; detector D2: Ax2 = 0.02 A, 
Wi = 1 X 10 -3 mc 2 ; the dotted line indicates the traversal time deduced from 
classical relativistic mechanics; (a) detector height W\ = 1 X 10~ 3 mc 2 ; (b) 
detector width Axi = 0.5 A 



shifted to higher traversal times. 

In the last part of this section, we fix the parameters of D\ at Ax\ = 0.5 A 
and Wi = 1 x 10~ 3 mc 2 and we vary the parameters Ax 2 and W2 of detector 
D 2 . We examine the following pairs of detector parameter Ax 2 = 0.02 A/W2 = 
1 x 10~ 3 mc 2 , Ax 2 = 0.02 k/W 2 = 1.0 mc 2 , Ax 2 = 0.5k/W 2 = 1 x 10" 3 mc 2 , and 
Ax 2 = O.5A/W2 = 1.0 mc 2 . We find that the resulting probability densities po and 
expectation values T tj o are nearly the same in the first three cases. The only exception 
is the case of a very wide and "height" detector (last case). In that case, the mean 
traversal time T tj o is lower than in the other cases. 

The probability Poo, 12 grows significantly if one increases the detector width Ax 2 
or the detector height W 2 . We get the same qualitative dependence of Poo, 12 on the 
parameters of detector D 2 as on the parameters of detector D\. 

Note that the following fact is true in the case of weakly intrusive detectors 
(Wi = W 2 = 1 x 10 -3 mc 2 ): the dependence of T ti0 on Ax\ is "stronger" than the 
dependence on Ax 2 . The reason for this is clear: changing the width Axi of the first 
detector D\ change not only the first "detection-time" but also the form of the wave 
function after the first detection. 

Summarizing, it exists a wide range of parameters of D\ and D 2 for which the 
mean traversal time does not change significantly. Remember that the same result 
was found in the study of time-of-arrival. 

6. Conclusion 

In we have proposed an extension of EEQT describing one spin i-particle in a 
relativistic, covariant way. 

In this paper, we have focused on applications of this formalism. We have 
calculated detection times of the particle in two-dimensional spacetime. The particle 
has moved freely except of the influence exerted on it by the detector (s). 

As a first application, we have computed the time at which the particle arrives at a 
detector ("free time-of-arrival"). We have found out that there exists good agreement 
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between the expectation values of our simulation and the results obtained by using 
classical relativistic mechanics of a free point particle. Moreover, we have shown that 
this fact is independent of the reference frame: the agreement is equally good in all 
reference frames. We have considered the situation with different detectors. We have 
found out that the good agreement between the results of our algorithm and the results 
obtained by using classical relativistic mechanics is not limited to a special value of 
detector parameters. It is obtained for a wide range of detector parameters. 

As a second application, we have examined the time difference between two 
detections of the particle by two detectors, one behind the other ("free traversal 
time"). Again, we have obtained good agreement between the simulated results and 
those obtained by classical relativistic mechanics in all reference frames. Moreover, 
we have shown that there again exist a wide range of detector parameters for which 
the mean traversal time does not change significantly. 

Summarizing, we have found interesting and meaningful results in these 
applications of our relativistic extension of EEQT. So we think that our formalism 
will be also useful in future applications, for example in the case, in which the particle 
is submitted to the action of a potential barrier. 
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